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We provide an analytic method to discriminate among different types of black holes on the ground 
of their strong field gravitational lensing properties. We expand the deflection angle of the photon in 
the neighbourhood of complete capture, defining a strong field limit, in opposition to the standard 
weak field limit. This expansion is worked out for a completely generic spherically symmetric 
spacetime, without any reference to the field equations and just assuming that the light ray follows 
the geodesies equation. We prove that the deflection angle always diverges logarithmically when 
the minimum impact parameter is reached. We apply this general formalism to Schwarzschild, 
Reissner-Nordstrom and Janis-Newman-Winicour black holes. We then compare the coefficients 
characterizing these metrics and find that different collapsed objects are characterized by different 
strong field limits. The strong field limit coefficients are directly connected to the observables, such 
as the position and the magnification of the relativistic images. As a concrete example, we consider 
the black hole at the centre of our galaxy and estimate the optical resolution needed to investigate 
its strong field behaviour through its relativistic images. 
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I. INTRODUCTION 

Gravitational lensing is one of the first applications of 
General Relativity ever studied 0]. Firstly it was recog- 
nized in the light deflection by sun, secondly in lensing 
of quasars by foreground galaxies, then in the formation 
of giant arcs in galaxy clusters and finally in galactic 
microlensing. Now it is an ordinary phenomenon in the 
panorama of astronomical observations (see || for a com- 
plete treatment and references therein). 

The full theory of gravitational lensing has been de- 
veloped following the scheme of the weak field approxi- 
mation and, in this formulation, it has been successfully 
employed to explain all the physical observations. 

In the last years, however, the scientific community 
is starting to look at this phenomenon from the oppo- 
site point of view, opening a strong field perspective. 
Viergutz |j| made a semi-analytical investigation about 
geodesies in Kerr geometry; in Ref. 0] the appearance 
of a black hole in front of a uniform background was 
studied; Falcke, Melia and Agol Q considered the emis- 
sion of the accretion flow as source. Virbhadra & Ellis 
H showed that a source behind a Schwarzschild black 
hole would produce one set of infinite relativistic im- 
ages on each side of the black hole. These images are 
produced when a light ray with small impact parameter 
winds one or several times around the black hole before 
emerging. Later on, by an alternative formulation of the 
problem, Frittelli, Kling & Newman |jj attained an exact 
lens equation, giving integral expressions for its solutions, 
and compared their results to those by Virbhadra & El- 
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lis. The same problem has been investigated by Bozza 
et al. in Ref. ||, where a strong field limit was first 
defined in Schwarzschild black hole lensing and used to 
find the position and the characteristics of all the rela- 
tivistic images analytically. Eiroa, Romero and Torres 
H applied the same technique to a Reissner - Nordstrom 
black hole. Recently, in another work pc| , Virbhadra 
& Ellis distinguished the main features of gravitational 
lensing by normal black holes and by naked singularities, 
analyzing the Janis, Newman, Winicour metric. They 
remarked the importance of these studies in providing a 
test for the cosmic censorship hypothesis. 

The reason for such an interest in gravitational lensing 
in strong fields is that by the properties of the relativis- 
tic images it may be possible to investigate the regions 
immediately outside of the event horizon. High resolu- 
tion imaging of black holes by VLBI could be able 
to detect the relativistic images and retrieve information 
about strong fields stored within these new observables. 
Moreover, since alternative theories of gravitation must 
agree with GR in the weak field limit, in order to show 
deviations from GR it is necessary to probe strong fields 
in some way. Indeed, deviation of light rays in strong 
fields is one of the most promising grounds where a the- 
ory of gravitation can be tested in its full form. 

Of course, the study of null geodesies in strong fields 
is not easy and up to now it has always been carried 
out using numerical techniques. An analytical treatment 
would enlighten the dependence of the observables on the 
parameters of the system, allow easy checks about the 
detectability of the images and open the way to compar- 
isons between the results in different metrics. In Ref. 
S, a new way to expand the deflection angle in the 
Schwarzschild metric was suggested. The deflection an- 
gle near its divergence was approximated by its leading 
order and its first regular term and then plugged into the 
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lens equation. In this way, very simple and reliable an- 
alytical formulae were derived for the relativistic images 
and their main features. 

So, as the weak field limit takes the first order devi- 
ation from Minkowski, the strong field limit starts from 
complete capture of the photon and takes the leading 
order in the divergence of the deflection angle. 

The strong field limit of Ref. || was only developed 
in Schwarzschild spacetime. In this paper, we provide 
a general method to extend the strong field limit to 
a generic static spherically symmetric spacetime. Our 
method is universal and can be applied to any spacetime 
in any theory of gravitation, provided that photons sat- 
isfy the standard geodesies equation. The parameters of 
the strong field limit expansion are directly connected 
with the observables, providing an effective tool to dis- 
criminate among different metrics. In Sect. 2, we state 
the problem and carry out the strong field limit of the de- 
flection angle. In Sect. 3, we apply the method to some 
simple metrics: Schwarzschild, Reissner-Nordstrom and 
Janis-Newman-Winicour black hole, discussing their dif- 
ferences with reference to the gravitational lensing phe- 
nomenology. In Sect. 4, we establish a connection be- 
tween the strong field limit coefficients and the relativis- 
tic images, analyzing the case of the black hole at the 
center of our galaxy as a concrete example where our 
results can be tested. Finally, Sect. 5 contains the sum- 
mary. 



II. STRONG FIELD EXPANSION OF THE 
DEFLECTION ANGLE 

A generic spherically symmetric spacetime has line- 
element [[l3| 

ds 2 = A{x)dt 2 - B{x)dx 2 - C{x) (dO 2 + sin 2 9<p 2 ) . (1) 
where 

A(xf=** 1 - 



B{xf^- 
C(xf 



1+ 2M 



(2) 



in order to correctly match the weak gravitational field 
far from the lensing object. 
We require that the equation 



C'{x) _ A'{x) 
C(x) ~ A(x) ' 



(3) 



admits at least one positive solution. We shall call the 
largest root of Eq. (Q) the radius of the photon sphere 
x m (for an alternative definition of the photon sphere, see 
Ul). A, B, C, A' and C must be positive for x > x m 

For metrics expressed in standard coordinates (C (x) = 
x 2 ) a sufficient condition for the existence of x m is the 
presence of a static limit (a radius x s such that A(x s ) — 



0) . Our strong field expansion takes the photon sphere as 
the starting point. In our study, therefore, we shall not 
consider naked singularities without a photon sphere. For 
a numerical study of this situation, see Ref. (lOj . 

A photon incoming from infinity with some impact pa- 
rameter u, will be deviated while approaching the black 
hole. It will reach a minimum distance xq and then 
emerge in another direction. Of course, the approach 
phase is identical to the departure phase, with the time 
reversed. 

By the conservation of the angular momentum, the 
closest approach distance is related to the impact pa- 
rameter u by 



(4) 



where the subscript indicates that the function is eval- 
uated at Xq. 

From the geodesies equation it is easy to extract the 
quantity 



B 



dx 



Co A 



1 



(5) 



which gives the angular shift of the photon as a function 
of the distance from the center (see [^3) for the complete 
derivation) . 

The deflection angle can then be calculated as a func- 
tion of the closest approach: 



a ( x o) = I ( x o) - 71" 

oo 

2\/Bdx 



c_Aa. _ i 

Co a 



(6) 
(7) 



It is easy to check that for a vanishing gravitational field 
(A = B = 1, C = x 2 ) a(xo) identically vanishes. In 
the weak field limit, the integrand is expanded to the 
first order in the gravitational potential. This limit is 
no longer valid when the closest approach distance sig- 
nificantly differs from the impact parameter (which, by 
Eq. (|j), means that A(xq) significantly differs from 1 
and/or C(xq) from Xq, i.e. the photon passes in a strong 
gravitational field). 

When we decrease the impact parameter (and con- 
sequently xo), the deflection angle increases. At some 
point, the deflection angle will exceed 2ir, resulting in a 
complete loop around the black hole. Decreasing u fur- 
ther, the photon will wind seve ral t imes before emerging. 
Finally, for xq = x m (see Sect. |ll A| ), corresponding to an 
impact parameter u — u m , the deflection angle diverges 
and the photon is captured (see Fig. |l|). 

In this section, we will show that this divergence is 
logarithmic for all spherically symmetric metrics. Our 
purpose is to get an analytical expansion of the deflection 
angle close to the divergence in the form 

a(x ) = -a log ( — - 1 ) + b + O (x Q - x m ) , (8) 
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FIG. 1: General behaviour of the deflection angle as a function 
of the closest approach xo. The deflection angle increases 
as xo decreases and diverges at xo = x m . Each time a(xo) 
reaches a multiple of 2ir, the photon completes a loop before 
emerging. 



where all the coefficients depend on the metric functions 
evaluated at x m . 

As the angular separation of the image from the lens 
is 9 = jj^- , where Do l is the distance between the lens 
and the observer, we need to express the deflection angle 



in terms of this variable. 
Eq.® to 



a{9) — —a log 



(QD 



OL 



1 + b + O (u - U m ) (9) 



which we define as the strong field limit of the deflection 
angle. The rest of this section is devoted to the calcu- 
lation of the two coefficients a and & in a first step and 
then of a and b. 



The function R(z, xo) is regular for all values of z and 
xo, while f(z,xo) diverges for z — > 0. To find out the 
order of divergence of the integrand, we expand the ar- 
gument of the square root in f(z,xo) to the second order 



/ (z,x ) ~ fo(z,x ) 



\/ az + (3z 2 



" =^r^(C' yo-C A' ) 







C A' 



(15) 
(16) 



(C C ' - 2C 2 ) y A' - C C yoA$] . (17) 



When a is non zero, the leading order of the divergence 
in /o is z -1 / 2 , which can be integrated to give a finite 
result. When a vanishes, the divergence is z~ x which 
makes the integral diverge. Examining the form of a, we 
see that it vanishes at xq = x m , with x m defined by Eq. 
(Q). Each photon having xo < x m is captured by the 
central object and cannot emerge back. 

To solve the integral (112) , we split it into two pieces 



I(x Q ) = I D (x Q ) + Ir{x ), 



So, we shall finally transform where 



Id(xo) = / R(0,x m )fo(z,x o )dz 



(18) 



(19) 



contains the divergence and 



Ir(xq) = / g(z,x )dz 



(20) 



A. Divergent term of the deflection angle 

We define two new variables 

y = A{x) 
y - yo 



l-2/o 



(10) 

(ii) 



where yo = A . The integral (0) in the deflection angle 
becomes 



I(x ) = / R(z,x a )f(z,x )dz 



R(z, x ) = 
f{z,x ) = 



CA!~ ^ ~ Vo ^ V °° 
1 

yyo- [{i-yo)z + yo]% 



(12) 

(13) 
(14) 



where all functions without the subscript are evaluated 
at x = A~ x [(1 -yo)z + y ]. 



g(z, x ) = R(z, x )f(z, x ) - R(0, x m )f (z, x ). (21) 

is the original integral with the divergence subtracted. 
We shall solve both integrals separately and then sum up 
their results to rebuild the deflection angle. Here we deal 
with Id and its divergence, while in the next subsection 
we shall verify that Ir is indeed regular. 

The integral Id(xq) can be solved exactly, giving 



I D {xo) = R{0, x m )—7= log -= . (22) 

VP V" 

Since we are interested in the terms up to O(xo — x m ), 
we expand a as 

a = 2 ^ mAm (x - x m ) + O(x - x m ) 2 , (23) 



l-y* 



where 

Pm — P\xQ—x ri 



C m (1 — ym) (C'^ym — C m A" (x m )) 



2y 2 nC m 



(24) 
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and substitute in Id(xq). Rearranging all terms, we find 



C. From a(xo) to a{9) 



Id(x ) = -a log 



xq 



R(0, Xj, 



, fl(0,i m ). 2(1 

b D = 7==— log 



1 ) + b D + O(x - x m ), (25) 

(26) 
(27) 



A m x m 



Eq. ( p5| ) yields the leading order in the divergence of 
the deflection angle, which is logarithmic, as anticipated 
before. The coefficient a of Eq. (|8|) is then given by Eq. 



B. Regular term of the deflection angle 

In order to find the correct coefficient b in Eq. (H) , 
we have to add to the term bo coming from Eq. (|27|), 
an analogous term coming from the regular part of the 
original integral, defined by Eq. j20|). 

We can expand I r (xq) in powers of (xq — x m ) 



Ir{xo) 



i 



n=0 



d n g 
dxVt 



dz (28) 



and evaluate the single coefficients. 

If we had not subtracted the singular part from 
R(z, XQ)f(z, xq), we would have an infinite coefficient for 
n = 0, while all other coefficients would be finite. How- 
ever, the function g(z, xq) is regular in z = 0, xq = x m 
as can be easily checked by a power expansion, recalling 
that a m — 0. 

Since we are interested to terms up to 0(xq — x m ), we 
will just retain the n = term 

l 

Ir(xo) = J g(z,x m )dz + O{x - x m ) (29) 
o 

and then 

b R = I R {x m ) (30) 

is the term we need to add to bo in order to get the 
regular coefficient. Recalling also the term — 7r in the 
deflection angle, we have 

b = -n + b D +b R . (31) 

The coefficient b R can be easily evaluated numerically 
for all metrics, since the integrand has no divergence. 
However, in many cases it is also possible to build a 
completely analytical formula for b R as well. In fact, in 
Schwarzsc hild m etric, the integral ( p9|) is solved exactly 
(see Sect. HIA). Then, in most metrics, we can expand 
Eq. (p9|) in powers of their parameters, starting from the 
Schwarzschild limit, and evaluate each coefficient sepa- 
rately. This is what we shall do for Reissner-Nordstrom 



From Eq. (g), we see that the minimum impact pa- 
rameter is 




(32) 



Expanding Eq. (|J), we find 

U - U rn = c(x - X m ) 

C" ii — C A" 




Using this relation, we can write the deflection angle 
as a function of 9: 



a{8) = —a log 



6D, 



OL 



u„ 



_ R(0,x m ) 

— lOg < —2±- ; 

2 Urn 



a lot 



20T, 

Urn 



(35) 
(36) 
(37) 



This concludes our general discussion of the form of the 
deflection angle in the strong field limit. Even if the proof 
is somewhat tricky, the application to concrete cases is 
very straightforward, as we shall see in Sect. [II. In fact, 
once we write a metric, it is sufficient to: 



and i?(0, x rn ) from Eq. 



and Janis- Newman- Winicour black holes (see Sect. Ill 



• solve Eq. (0) to find x n 

• write [3 m from Eq. (| 

©; 

• compute b R from Eq. (|3(]) numerically or by a 
proper expansion in the parameters of the metric; 

• compute the coefficients u m , a and b from Eqs. 
©, (If) and @, respectively. 

The crucial step is the calculation of b R , since it is the 
only integral involved in the whole procedure. 



III. APPLICATIONS 

In this section, we apply the general formulation of 
the strong field limit to three simple examples. First, we 
shall revisit Schwarzschild spacetime, whose strong field 
limit gravitational lensing has already been studied in 
Ref. B. Then we will apply our method to Reissner- 
Nordstrom spacetime, which was explored numerically 
in Ref. || . Finally, as an example of extended theory of 
gravitation, we shall consider a Janis, Newman, Winicour 
black hole, whose strong field limit expansion has not 
been investigated so far. In all three cases, we shall derive 
analytical formulae for the strong field limit coefficients, 
in order to analyze the functional dependences of the 
deflection angle on the parameters of the metrics. 
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A. Schwarzschild lensing 

This is the simplest spherically symmetric metric de- 
scribing the outer solution for a black hole. It only de- 
pends on the mass of the central object (by Birkhoff's 
theorem). It is convenient to define the Schwarzschild 
radius x s = 2M as the measure of distances; then, in 
standard coordinates, the functions in the metric become 



A(x) = 1 - 



1 



B(x) = ( 1 - 
C(x) = 



(38) 

(39) 
(40) 



which obviously satisfy all hypotheses required in Sect, 
[ij, with static limit x s = 1. 

The two functions R(z,xq) and f(z,Xo) read 



R(z,x ) = 2 
f(z,x ) = 



\( 2 ') 



(41) 
,(42) 



FromEqs. ©, @, or directly from the expansion of 
the denominator of /, we read off the coefficients a and 

P 



a 



2-1 

Xq 

i. 

Xq 



(43) 
(44) 



The equation a = defines the radius of the photon 
sphere 



(45) 



Consequently, (3 m = 1. 

In this simple case, it is possible to solve the integral 
in Eq. (|3(]) exactly and write the regular term in the 
deflection angle 



6 fi = 21og 6(2-Vtj 



= 0.9496. 



(46) 



From Eqs. ( pq ) and (p7[), we derive the coefficients a, 
6 and u m of the deflection angle 



a= 1 

b = — 7T + 6_R 

3^ 



log 6 = -0.4002 



(47) 
(48) 

(49) 



Then the Schwarzschild deflection angle, in the strong 
field limit, is 



a(0) = -log 



26D 



Oh 



3\/3 



log 



216 (7-4V3 



— 7T, 

(50) 



as obtained in Ref. [||. 

In this simple case, we can compare the exact deflection 
angle a ex (8), calculated numerically and the strong field 
limit a app (9). The most external image appears where 
a(0) falls below 2tt. This happens at u — u m = 0.003264. 
Here the discrepancy between a ex and a app is about 
0.06%, which corresponds to an error in the position of 
the outer image of the order 0.4%. With such a high 
accuracy, we are encouraged to take the Schwarzschild 
strong field limit as the starting point for successive se- 
ries expansions to evaluate bn in more advanced metrics. 



B. Reissner-Nordstrom lensing 

Reissner-Nordstrom metric describes the gravitational 
field of a spherically symmetric massive object endowed 
with an electric charge q. The metric functions in stan- 
dard coordinates are 



A(x) 

B(x) 
C(x) 



1 



1 

x 



1 

x 



(51) 

(52) 
(53) 



They satisfy the hypotheses required in Sect. [n|, only 
when q < However, beyond the critical value q — 

0.5, there is no event horizon and causality violations 
appear @ H. We shall restrict to q < 0.5. 
The coefficients a and (3 are 



3 

x 



x 



xl J x - 2q' 2 



\xo 4 xl) (xl-2q>f 



(54) 
(55) 



which reduce to the Schwarzschild coefficients when q — > 
0. From the equation a = 0, we derive the radius of the 
photon sphere 



T U + \ 1 



32q 2 



(56) 



which yields 



9 + 32g 2 - 144g 4 + b\2q & + ^9 - 32q 2 
(3 + 16<? 2 - 80g 4 )] [8 (q - 4q 3 )] ~ 2 . (57) 



The regular term bn cannot be calculated analytically, 
but we can expand the integrand in ( p9| ) in powers of q 
and evaluate the single coefficients. We get 



b R = b R , a + b Ra q 2 + 0(q 4 ), 



(58) 



where b R is the value of the coefficient for an uncharged 
black hole, calculated in the previous subsection and 
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given by Eq. (flq); the correction is quadratic in the 
charge q of the black hole. Its coefficient is 



AR,2 



{x/3 



log 



-1.5939. 



(59) 

It is very easy to calculate further terms in the expan- 
sion of bR, deriving analytical formulae which prove to 
be very accurate even for large values of q. 

Following Sect. ||, we calculate the coefficients for the 
formula of the deflection angle 



- l in 



\/(3 - x m ) x 2 n - 9q 2 x m + 8 

= — 7r + bn + a log 2 (x m - q 2 ) 

[(3 - x m ) x 2 m - 9q 2 x m + 8g 4 ] 

(x m - 2q 2 f (x 2 m - x m + q 2 ) 
2 



3 + y/9 - 32q< 



4V2a/3 - 8q 2 + V9 - 32q 2 



(60) 

(61) 
(62) 



In H , the coefficients a and b were calculated numeri- 
cally. Here, by our general method, we have been able to 
derive the coefficient a for the logarithmic divergence ex- 
actly and find a formula for b which is valid up to second 
order in q, indicating the way to extend it to an arbitrary 
order. 

We notice that the radius of the photon sphere x m 
decreases as the charge increases, but becomes imaginary 
only for q — > i; i.e., even when there is no horizon, 

a photon can be captured by the gravitational field of a 
hypothetical object with a charge greater than the critical 
value. 

In Fig. |^, we evaluate the deflection angle at u — 
u m + 0.003 as a function of q. The plot shows that, 
just using the first correction to bn, we get an excellent 
approximation to the deflection angle. In fact, we see 
that up to q = 0.3, the error in the position of the outer 
image, calculated using the strong field limit, stays below 
4%. 

Finally, in Fig. ||, we plot the coefficients of the strong 
field limit as functions of q, calculating the full b numeri- 
cally. We see that a and b deviate from the corresponding 
Schwarzschild values as the charge increases. As we shall 
see in Sect. IV, the strong field limit coefficients are di- 



rectly connected to the observables. It is then possible, 
in principle, to distinguish a Reissner-Nordstrom black 
hole from a Schwarzschild black hole, using strong field 
gravitational lensing. 



C. Janis-Newman-Winicour lensing 



The spherically symmetric solution to the Einstein 
massless scalar equations (i? M „ = $:^) can be 

written in Janis-Newman-Winicour (JNW) coordinates 
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FIG. 2: Deflection angles in Reissner-Nordstrom metric eval- 
uated at u — u m + 0.003 as functions of q. The solid line is 
the exact deflection angle; the dashed line is the strong field 
limit with bn truncated to second order in q. 
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FIG. 3: Coefficients of the strong field limit in Reissner- 
Nordstrom metric as functions of q. 



H 



A(x) 
B{x) 
C{x) 



= 1- 



1 



7 



X 
1 

1 

x 



1-7 



2y/AP + q 

M 



11 f 10 ^ 1 " 1 

2 \ X 



7 : 



y/M 2 + q 2 



(63) 
(64) 
(65) 
(66) 
(67) 



wh ere all d istances are measured in terms of x s = 
2\J M 2 + q 2 and q is the scalar charge of the central ob- 
ject. This metric admits a photon sphere external to the 
static limit when 7 > 4, i.e. when q < M. We shall thus 
restrict our investigations to objects with scalar charge 
lower than their mass. In Ref. 1 10 , the gravitational lens- 
ing of this object was investigated numerically even when 
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FIG. 4: Deflection angles in a JNW black hole evaluated at 
u = M m + 0.003 as functions of 7. The solid line is the exact 
deflection angle; the dashed line is the strong field limit with 
bn truncated to first order in 7. 



q > M . In this situation, it was shown that a drastically 
different and interesting phenomenology shows up. 

As in the previous cases, we compute all the coeffi- 
cients, taking into account that our metric is not written 
in standard coordinates. The coefficients a and (3 are 



V xo ) lx l 
(2 7 + 1) (7 + 1) - 2a;o (3 7 + 1) + 



(x - 1)" 



2^x1 



(68) 



(69) 



which easily reduce to Eqs. ( [4 3D and ( ]44| ) when 7=1. 
From the Eq. a = 0, we derive the radius of the photon 
sphere 



and then 



2 7 +l 



[(2 7 +i) 7 -(2 7 -ir 

4 7 2 ( 4 7 2 _ ^27-1 



(70) 



(71) 



In the same way as in Reissner-Nordstrom metric, in- 
creasing the charge of the central object, the photons 
are allowed to get closer to the centre. When 7 = 0.5, 
x m = %s and we should change our coordinates frame to 
perform any significant study. 

Here as well we cannot solve the integral (p9| ) exactly. 
However, we can expand the integrand in powers of (7— 1) 
to get 



b R = b° R - 0.1199(7 - 1) + 0(7 - l) 2 



(72) 



FIG. 5: Coefficients of the strong field limit in a JNW black 
hole as functions of 7. 



Finally we compute the coefficients of the strong field 
limit 



a = 1 

b = — 7T 



+2 log 

(2 7 



-b R + 

[(2 7 + l) 7 -(2 7 -l) 7 ] 2 (27 + l) 



2 7 2 (2 7 - 1) 



27-1 



2(27-1) 



7-5 



(73) 

(74) 
(75) 



with b R given by Eq. (p 



Surprisingly, the leading order coefficient a is the same 
as in the Schwarzschild case, in spite of the fact that the 
radius of the photon sphere has changed. 

Once we fix u — u m = 0.003, we see in Fig. |i| that 
the deflection angle decreases as we increase the charge 
(decrease 7) until 7 reaches the value 0.6, then the deflec- 
tion angle increases again. The strong field limit, with 
bn truncated to first order in 7, at 7 = 0.7 is precise up 
to 4% in the determination of the outer image. 

In Fig. ^, we plot the coefficients of the strong field 
limit, a is constant, b increases as 7 decreases, but u m 
decreases enough to make the deflection angle decrease 
at constant u — u m until 7 reaches the value 0.6. Com- 
paring with Fig. ||, we can observe that a JNW charge 
has completely different effects on the strong field limit 
coefficients than an electric charge and can be identified 
without confusion. 



IV. OBERVABLES IN THE STRONG FIELD 
LIMIT 

In Sect. H we have proved that the strong field limit 
approximation can be used to obtain a simple and reli- 
able formula for the deflection angle, which contains a 
logarithmic and a constant term. Now we plug the for- 
mula (||) into the lens equation and establish direct rela- 
tions between the position and the magnification of the 
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relativistic images and the deflection angle, calculated 
according to the strong field limit. 

The lens equation in the strong field limit was derived 
in Ref. il. It reads 



Dls 
Dos 



(76) 



where Dls is the distance between the lens and the 
source, Dos — Dol + Dls, P is the angular separation 
between the source and the lens, 9 is the angular separa- 
tion between the lens and the image, Aa n = a{9) — 2mr 
is the offset of the deflection angle, once we subtract all 
the loops done by the photon. 

To pass from the deflection angle a(9) to the offset 
Aa n , we need to find the values 0° such that a(9^ l ) = 
2mt. Solving Eq. (35) with a(9) = 2mr, we find 



Dol 



(l + e n ) 



(77) 
(78) 



The offset Aa n is then found expanding a{9) around 
e = 9 a n . Letting A9 n = 9 - 0° we find 



a aD L AQ 
Aa n = A9 n . 

The lens equation becomes 



(79) 



/ 3=(0°+A0„) + (^^W (80) 

V u rn^n D S J 

Now we derive the position of all relativistic images, 
their magnification and the critical curves of the lens. 

The second term in the r.h.s. of Eq. ( |S0| ) is negligi- 
ble when compared to the last one (since u m *C Dol)- 
Immediately, we find 



u m e n ((3 - 91) Dps 
aD L sD L 



(81) 



where the correction to 0° is much smaller than 0„ . 

This formula is valid both for the images on the same 
side of the source and the images on the opposite side. In 
fact, to find the latter, it is sufficient to take a negative 
13 in Eq. @. 

Finally, we have expressed the position of the relativis- 
tic images in terms of the coefficients a, b and u m . If we 
manage to determine these coefficients from the observa- 
tion of the relativistic images, we are rewarded with in- 
formation about the parameters of the black hole stored 
in them. 

The critical curves are just Einstein rings correspond- 
ing to a source perfectly aligned with the lens. Their 
radius is obtained putting j3 = in (|3l|). 

Another important source of information is the magni- 
fication of the images which is the inverse of the Jacobian 
evaluated at the position of the image. For simplicity, we 



approximate the position of the images by 0„, recalling 
that the correction provided by Eq. (J8l|) is negligible. 



1 



/<-« 



pop 
o do 



We have 



dp 
89 



= 1 



aDpL Dls 
u m e n Dos 



(82) 



(83) 



where the first term is small compared to the second and 
can be neglected. 

Finally, the magnification is 



1 



0° 



\detJ\go 



<>l> 



u m (1 + e n ) Dos 
1 aPD 2 OL D LS 



, (84) 



which decreases very quickly in n. 

Formulae (|77]) and (|84|) relate the position and the 
magnification to the strong field limit coefficients. The 
successive step is to solve the inversion problem, i.e. we 
have to find out the most effective way to go back from 
measured positions and fluxes to the strong field limit co- 
efficients, which carry the information about the nature 
of the black hole. 

The minimum impact parameter can be simply ob- 
tained as 



u m = Dol9 



(85) 



where (9oo represents the asymptotic position approached 
by a set of images, obtained by Eq. ( |77| ) in the limit 
n — > oo. 

To obtain the coefficients a and 6, we need to separate 
at least the outermost image from all the others. We 
shall thus consider the simplest situation where only the 
outermost image 9\ is resolved as a single image, while 
all the remaining ones are packed together at #oo. 

Our observables will thus be 



s = 9 1 



n=2 



(86) 
(87) 



which respectively represent the separation between the 
first image and the others, and the ratio between the flux 
of the first image and the flux coming from all the other 
images. 

The sum of the fluxes of all the set of relativistic images 
except the first is 



u; n D se- e- 



a(3D OL D LS - 1 



(88) 



We can notice that 3> 1 and e« is of order one, 
since a and b are of order one too. Using these simple 
observations, we can simplify our formulae to have 
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a = O^ei-^ (89) 
r = e*. (90) 

These two formulae can be easily inverted to give 



6 = alog^— y (92) 

Finally, just measuring an angular separation and a flux 
ratio, we are able to reconstruct the full strong field limit 
expansion of the deflection angle for the observed gravi- 
tational lens. 

The coefficients of the strong field limit are constrained 
by the characteristics of the metric to be well precise 
real numbers. Thus, if strong field gravitational lensing 
will be detected, comparing the experimental coefficients 
with the theoretical expectations, calculated according 
to different models, we are able to identify the nature of 
the lensing black hole unambiguously. 

A. A numerical example: lensing by the galactic 
supermassive black hole 

It is significant to consider a realistic case where we 
can discuss the instrumental sensitivity required to de- 
tect relativistic images and possibly distinguish between 
different black holes through the reconstruction of the 
strong field limit coefficients. 

The centre of our galaxy is believed to host a black 
hole with mass M = 2.8 x 10 6 M Q Q. The lensing 
of a background source by this supermassive black hole 
was discussed in detail by Virbhadra & Ellis || . Taking 
Dol = 8.5kpc as the distance between the sun and the 
center of Galaxy, they found that the separation between 
each set of relativistic images with respect to the central 
lens would be #oo ~ 17 microarcsecs. In principle, such 
a resolution is reachable by actual VLBI projects, but 
we must be aware that the disturbances intrinsic in such 
observations (mainly due to extinction and emission by 
accreting matter), would make the identification of the 
relativistic images very difficult, as already pointed out 
inRef. §. 

In Tab. following this line, we estimate the quantities 
we need for a complete strong field limit reconstruction in 
different situations, starting from a simple Schwarzschild 
black hole and then going to black holes with different 
values of the electric charge and JNW charge. 

Looking at Tab. |, we can make some considera- 
tions of different order. Indeed, the easiest parame- 
ter to evaluate is the minimum impact parameter u m , 
since a microarcsec resolution is reachable in the next 
years. This information alone, can already distinguish 
between a Schwarzschild or other types of geometry. In 



fact, since the total mass and the distance of the black 
hole are known to a reasonable accuracy (and possi- 
bly will be even better fixed in the next years), an VLjn 
smaller than predicted would signal that the structure 
of spacetime close to the central black hole is not de- 
scribed by Schwarzschild solution. On the other hand, if 
u m is compatible within experimental uncertainties with 
Schwarzschild case, we would set an upper limit for the 
parameters describing other black holes, such as an elec- 
tric or a scalar charge. 

In a second extent, to fit all the strong field limit co- 
efficients into any black hole model, we need to separate 
at least the outermost relativistic image from the oth- 
ers. We see that this can be done only increasing optical 
resolution at least by two orders of magnitude with re- 
spect to actual observational projects. Therefore, with 
these numbers, it seems that we are forced to wait for 
further technological developments. However, given the 
evolution rate of astronomical facilities in the last twenty 
years, it is not unthinkable that these two orders of mag- 
nitude will be reached within a not so far future. 

Black holes are also present in the bulge of other galax- 
ies. As far as we have investigated, in the best cases, 
the instrumental resolutions needed for strong field grav- 
itational lensing are about the same as for the Milky 
Way black hole. So, the determination of u m by next 
future observations would possibly confirm or disprove 
Schwarzschild geometry for several extra-galactic black 
holes too. The complete determination of the strong field 
limit coefficients instead remains a long term project, un- 
less more favourable astrophysical situations emerge. 



V. SUMMARY 

Gravitational lensing is undoubtedly a potential pow- 
erful tool for the investigation of strong fields. By general 
arguments we have shown that the deflection angle di- 
verges logarithmically as we approach the photon sphere. 
We have drawn a general method to compute the coef- 
ficient of the leading order divergent term and the first 
regular term. When the latter cannot be calculated ana- 
lytically, we have seen that it can be well approximated 
by a simple series expansion starting from Schwarzschild 
spacetime. 

We have applied our method to Schwarzschild, 
Reissner-Nordstrom and Janis-Newman-Winicour black 
hole, explicitly calculating and plotting the strong field 
limit coefficients. 

Of course, it is possible to apply the strong field limit, 
in the form given in this paper, to any spherically sym- 
metric metric representing a black hole. In this way, it is 
possible to compare the gravitational lensing behaviour 
of these objects in different theories of gravitation. In 
principle, the extension to non-sphcrically symmetric and 
rotating black holes is possible. However, the dependence 
of the deflection angle on more than one variable can put 
severe obstacles in the way of analytic solutions of the 
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Schwarzschild 


Reissner - 


Nordstrom 


Janis-Newman-Winicour 


q = 0.1 


q = 0.2 


q = 0.3 


q = 0.4 


7 = 0.9 


7 = 0.8 


7 = 0.7 


7 = 0.6 


Boo (/xarcsecs) 


16.87 


16.76 


16.41 


15.78 


14.76 


16.67 


16.38 


15.93 


15.13 


s (/Ltarcsecs) 


0.0211 


0.0216 


0.0234 


0.0275 


0.038 


0.0213 


0.0216 


0.0222 


0.0239 


r m (magnitudes) 


6.82 


6.79 


6.69 


6.49 


6.07 


6.82 


6.82 


6.82 


6.82 


Um/Rs 


2.6 


2.58 


2.53 


2.43 


2.27 


2.57 


2.52 


2.45 


2.33 


a 


1 


1.005 


1.02 


1.052 


1.123 


1 


1 


1 


1 


b 


-0.4002 


-0.3993 


-0.3972 


-0.3965 


-0.4136 


-0.3808 


-0.35 


-0.2945 


-0.1659 



TABLE I: Estimates for the main observables and the strong field limit coefficients for the black hole at the center of our 
galaxy in different hypotheses for the spacetime geometry. 9^ and s are defined in the text; r m is r converted to magnitudes: 
r m = 2.5 Log r; u m , a and b are the strong field limit coefficients; Rs = 2l ^j v/ is the Schwarzschild radius. 



problem. Nevertheless, this is indeed another important 
point which needs to be investigated to complete the pic- 
ture of black hole lensing. 

Differences in the deflection angle are immediately re- 
flected on the relativistic images. If the mass and the 
distance of the lens is known, then the detection of any 
set of relativistic images would immediately check the 
Schwarzschild geometry. VLBI should be able to provide 
an observational answer, if the relativistic images are not 
hidden behind environmental noise. 

Furthermore, if the outermost image is resolved from 
the others, it is then possible to fully reconstruct the 
strong field limit coefficients and select a precise black 
hole model. Our present observational facilities are not 
so far from the required resolutions, which, for the galac- 
tic black hole, are of the order of 0.01 microarcsecs. As 
a long term project the detection of the outermost im- 
age stands as a very interesting, non-trivial challenge for 
future technology. 



Strong field limit represents an important step in the 
construction of a robust theoretical scheme connecting 
the gravitational lensing with the strong field properties. 
By a simple and reliable expansion, it clarifies the whole 
phenomenology and the differences between various mod- 
els that should be expected in the appearance of relativis- 
tic images. If these so elusive features will be detected, 
we will finally have a way to effectively discriminate be- 
tween alternative theories of gravitation and grow in our 
knowledge of spacetime. 
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